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Abstract. The two-nucleon interaction derived in chiral effective field theory is used to investigate the 
saturation properties of symmetric nuclear matter at next-to-leading order. Pauli blocking effects in the 
two-pion exchange diagrams are shown to produce an effective attraction in nuclear matter. 

PACS. 21.65.-f Nuclear matter - 21.65.Mn Equation of state of nuclear matter - 21.30.-x Forces in 
hadronic systems 



> 

00 

cn 

(N 
O 
OO 

O 



Within the last two decades, effective field theory has 
provided new methods for nuclear structure investigations. 
Two-, three- and four-nucleon interactions based on the 
most general chiral effective pion-nucleon Lagrangian have 
been developed based on Weinberg's power counting scheme. 
The phase shifts of the nucleon-nucleon interaction below 
the pion production threshold are reproduced with a preci- 
sion that is comparable with the one achieved by modern 
phenomenological two-nucleon potentials, when pushing 
the expansion to the order N'^LOQ see Refs.[l,2 . These 
forces together with three-nucleon forces at N'^LO have 
been applied to a cornucopia of reactions. In particular, 
various three-nucleon scattering observables have been an- 
alyzed at low and intermediate energies showing, in gen- 
eral, a good agreement with the data. The binding energies 
of ^He and ''He have been reproduced successfully once the 
leading three-nucleon force is included. Within the frame- 
work of the no-core shell model, the chiral forces have 
been used to study nuclei with A ~ 7, 10 ... 13 systems, 
in particular the sensitivity to the three-nucleon force, 
see [3111] for more details. Recent progress in these fields 
is reviewed in Ref. [5]. Given these successes, the ques- 
tion arises whether effective field theory can contribute to 
the solution of the nuclear matter problem. First ground- 
breaking attempts in this direction have been done by the 
Munich group, see e.g. Refs. (and references therein). 
In the present communication, we apply the NLO version 
of the chiral nucleon-nucleon interaction to symmetric nu- 
clear matter. This differs from the Munich approach which 
concentrates on the development of a chiral effective inter- 
action to be used exclusively in the nuclear medium with- 
out link to the two-nucleon interaction in the vacuum. 
In Ref. [8], nuclear matter has been studied with a low- 
momentum two-nucleon force Viow-k derived from the Ar- 



^ Throughout, we use the abbreviation LO for leading order, 
NLO for next-to-leading order and so on. 



gonne wis potential via renormalization group techniques. 
It was found that the low-momentum interactions can be 
used in second order perturbation theory to describe nu- 
clear matter, provided Pauli blocking effects due to the 
medium are properly included. This finding suggests that 
a resummation of the ladder diagrams, as it is mandatory 
in Brueckner theory, can be avoided. Recently, the simpli- 
fied version AV6' of the Argonne potential has been used 
in a quantum Monte Carlo study of symmetric nuclear 
matter [9 . For Fermi momenta less than kp = 1.33 fm~^, 
the binding energies obtained for nuclear matter are not 
smaller than the ones calculated within the Fermi hy- 
pernetted chain approach or the Brueckner Hartree-Fock 
method [10]. This fact is interpreted as evidence for the 
inadequacy of three-nucleon interaction models for all den- 
sities 12] ■ 

Both the two-nucleon potential derived from effective 
field theory and T4ow-k derived from any of the high- 
precision NN potentials can be characterized by the fact 
that the large momentum components of the interaction 
are removed, in contrast to meson-exchange potentials. 

The chiral potential differs from phenomenological ones, 
such as wi8, by its power counting scheme which allows sys- 
tematic improvements. In Ref.}13j. the pion-nucleon ver- 
tices are systematically taken from the chiral expansion 
of the pion-nucleon interaction. The two-nucleon inter- 
action is generated by eliminating the pion-two-nucleon 
states via the FST-Okubo projection method [TT1IT2] . The 
two-pion exchange diagrams relevant for the discussion of 
nuclear matter properties are shown in FigU] for the or- 
der NLO. The Weinberg- Tomozawa term generates a class 
of triangle diagrams (c) characterized by one internal nu- 
cleon line. Moreover, there are diagrams with two internal 
nucleon lines, such as the crossed-box two-pion exchange 
graph (a). Notice that the class of box two-pion exchange 
diagrams (b) in Fig[T] is not two-nucleon reducible and 
therefore has to be iterated when solving the Schrodinger 
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a) 



b) 



c) 



d) 



Fig. 1. Two-pion exchange corrections to the chiral two- 
nucleon potential in NLO. a) crossed box graph, b) wave func- 
tion renormalization diagram, c) Triangle graph, and d) foot- 
ball diagram. Solid/dashed lines denote nucleons/pions. The 
thick lines indicate medium modifications due to Pauli block- 
ing. 
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Fig. 2. The binding energy per particle of nuclear matter as 
a function of the Fermi momentum kp- Solid (red) line: The 
Pauli blocking effects due to Eqs. (|2l3p are included. Dashed 
(blue) line: Pauli blocking effects of Fig. 1 not included. Dash- 
dotted (violet) line: AFDMC results [9]. The empirical satura- 
tion point is given by the black circle. 



equation for two-nucleon scattering [14] . Also present are 
the so-called football diagrams (d), however, these do not 
have internal in- medium nucleon propagators. In the nu- 
clear medium, the two-nucleon interaction evaluated in 
the vacuum is strongly modified by Pauli blocking effects. 
The most important and well-known Pauli blocking effect 
affects the intermediate two-nucleon states. For the chiral 
two-nucleon interaction, there are additional Pauli block- 
ing effects due to the intermediate one- and two-nucleon 
propagators in the diagrams shown in FigU] Formally, 
these Pauli blocking effects can easily be taken into ac- 
count. The chiral potential is derived in Ref. [14 using the 
FST-Okubo method of unitary transformations, introduc- 
ing a projector 77 on the two-nucleon space and a projector 
A for the remaining part of the Fock space which consists 
of two-nucleon and one-pion or multi-pion states. The cor- 
responding unitary operator is parametrized in terms of 
the transition operator A = XAq which fulfills the decou- 
pling equation 

\{H -[A,H]- AHA)ri^Q. (1) 

This equation is solved perturbatively utilizing the chi- 
ral power counting. Once A is computed at a given order 
in the chiral expansion, the effective Hamilton operator 
acting on the nucleonic subspace of the Fock space can be 
obtained in a straightforward way. In the nuclear medium, 
one has to introduce medium-modified projectors 

fi = j d^h J d'fc2(i-e(fei))(i-e(fc2)) 

X I fcifc2)(fcifc2 I, (2) 

A = J d^'k, J d^'h J d^'q^il - 0(fci))(l - 0(fc2)) 

X I kik2qn){kik2qn | + . . . , (3) 

ellipses refer to components with two and more pions, 
0{k) = 1 for /c < fcF and 0{k) = for fc > kp- The 



momenta of the two nucleons and the pion are denoted by 
fci,fc2 and QTr, repectively. The diagrams shown in Fig[T] 
require loop integrations over non-trivial integration vol- 
umes defined by the Pauli blocking constraints. The bind- 
ing energy is obtained by summing all two-nucleon re- 
ducible diagrams, as it is done in Brueckner-Hartree-Fock 
theory. Therefore, the perturbative calculation of Ref. [8] 
can be considered as the limit of the present approach 
for weak potentials. At NLO, genuine three-body forces 
are not yet included0 The two-nucleon interactions de- 
rived in Ref. [2] are characterized by two regulators, the 
cut-off A introduced in the spectral function representa- 
tion of the two-pion exchange diagrams and the cut-off A 
required for the non-perturbative renormalization of the 
Lippmann-Schwinger equation. The actual potentials are 
derived for a large set of regulators A and A in order to 
generate theoretical error bands for a given order of the ex- 
pansion. Nuclear matter introduces a new scale into effec- 
tive field theory, the Fermi momentum, empirically given 
hy kp — 262 MeV. The effective two-nucleon interaction 
in nuclear matter should be able to describe momentum 
transfers of about 2k p. Therefore we have restricted the 
nuclear matter calculations to those potentials which have 
Lippman-Schwinger regulators A > 500 MeV. Fig.[2]shows 
the binding energy per particle as a function of the Fermi 
momentum kp for the NLO potential of Ref. [2^ derived 
for A = 550 MeV and A = 500 MeV. densities. One finds 
that at NLO, the chiral two-nucleon interaction produces 
a saturation point for nuclear matter, though not at the 
empirical values for the density and binding energy. 

In Refs. [H], the three-body forces were found to be 
essential in order to obtain saturation in a perturbative 
scheme. In the chiral approach of Ref. Pfj, the effect of 
three-body forces has been studied in detail. The calcu- 

A rough estimate of the effect of three-nucleon forces has 
been given in Ref. [16j . 
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lation including three-body forces saturates at the em- 
pirical density of nuclear matter p^^p = 0.160fm~^, if a 
pion-nucleon form factor of monopole type is used with 
a cut off Aq = 1145MeV. A calculation without three- 
body forces can obtain saturation at pexpi provided the 
cut off is reduced to Aq = 770 McV. The saturation found 
in the present calculation is entirely due to the regulators 
because there are no three-body forces at NLO. At order 
N^LO, chiral three-body forces appear which will reduce 
the sensitivity of the results to the choice of the regulators. 

The present calculation provides quantitative results 
concerning the relevance of the Pauli corrections imposed 
by Eqs. (|2|3p . We parameterize the saturation curves fol- 
lowing Ref. 9] by 

^^^+a{x^xf+P{x~xf, (4) 

with X = p/pcxp and pexp — O.lGOfm^'^ and obtain a 
binding energy Et^/A — — 18.35MeV at a Fermi momen- 
tum kp — 1.60 fm^ for the calculation which includes the 
Pauli blocking effects due to Eqs. (|2l3p . The compression 
modulus of nuclear matter then is given hy K — 279 MeV, 
corresponding to a = 5.0 MeV and (3 = -0.34 MeV. This 
is slightly larger than the value of the compression mod- 
ules derived from giant resonances, K = (260 ± 10) MeV 
[17] . A discussion of the model dependence of the compres- 
sion modulus can be found in Ref. TR' . If one discards the 
Pauli blocking effects due to Eq.Q, the saturation point 
moves to Eq/A = —15.30 MeV and kp — 1.57 fm^ , and a 
compression modulus oi K = 244 MeV is obtained. In or- 
der to understand the medium effects due to Eqs. (|2l3p . we 
have studied the influence of the various classes of the bare 
two-nucleon interactions on the two-nucleon phase shifts 
by multiplying the contributions of individual diagrams 
with some numerical factor. For the most important '^Si 
and ^Sq phase shifts, the (numerically) dominant contri- 
bution arises from diagrams a) in Fig[l] and is repulsive 
in both cases. In nuclear matter, the Pauli blocking of 
the diagrams involving two pion lines effectively reduces 
the repulsive contribution to the two-nucleon potential. 
Therefore, a gain of binding energy results, as seen in the 
calculations shown in Fig. [2l Quantitatively, the gain of 
binding energy is large enough to pull the saturation curve 
below the empirical saturation point of nuclear matter. 
The overall gain of binding energy in the vicinity of the 
saturation density is about 3 MeV. Genuine three-body 
interactions tend to yield an overall repulsive contribu- 
tion to the binding energy. The new class of Pauli block- 
ing effects studied in the present work may therefore be 
helpful to move the theoretical saturation point closer to 
the empirical value. Ref. ^ shows error bands for each 
phase shift which reflect the theoretical uncertainty due 
to the choice of the regulators. In Fig. [3] we show the cor- 
responding uncertainties for the nuclear matter saturation 
curve limiting ourselves to the variation of the cut-off A. 
The calculations have been performed omitting the Pauli 
blocking effects due to Eqs. (I2I3|) because the full inclu- 
sion of those effects is numerically time-consuming. This 
simplified calculation, however, already provides a reason- 
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Fig. 3. Dependence on the spectral function cut-off of the 
saturation curve. The Pauli blocking effects due to Eqs. ([2]) 
are included. The solid (black) , dashed (blue) and dash-dotted 
(green ) line refer to the spectral function cut-offs A — 500, 600 
and 700 MeV, respectively. 



able estimation of the theoretical uncertainty. In the limit 
of low densities one observes that the binding energies do 
not depend on the regulator employed. A similar obser- 
vation can be found in Ref. 7 . It may be understood as 
an illustration of Hugenholtz's well-known justification of 
Brueckner's method to describe the energy of a Fermi gas 
in the low density limit by summing the ladder diagrams 
of the two-body interactions 15 . At NLO, the saturation 
points are clearly cut-off dependent. The Fermi momen- 
tum at the saturation point varies from kp = 1.59 fm~'^ 
to kp = 1.66 fm~^, while the binding energy per particle 
moves from -15.3 MeV to -17.8 MeV. 

One of the most successful effective interactions to be 
used for the description of nuclear ground states has been 
proposed by Skyrme, see Ref. 18] for a recent review. In 
its simplest version, this interaction introduces contact in- 
teractions which have an analytical structure similar to 
the contact terms of the chiral two-nucleon potential. The 
Skyrme interaction differs from chiral interactions by in- 
troducing a contact three-body term 

w(ri,r2,r3) = t^Siri - r2)S{r2 - ra). (5) 

A term like this is not produced by the leading chiral 
three-body interactions jT6] . Moreover, if interpreted as a 
genuine three-body interaction, it leads to well-known spin 
instabilities ^19j in nuclei signalled by the Landau parame- 
ter go = —t^/S. In practical calculations, one therefore re- 
places the three-body interaction by a density-dependent 
two-body interaction v = t^/Gp"' with a = 1, effectively 
freezing the unphysical spin dependence of Skyrme's three- 
nucleon interaction. For the following comparison, we do 
not try to interprete the density-dependent term as a gen- 
uine three-body interaction, but rather as an effective rep- 
resentation of a regulator. In Ref. |14| . an exponential reg- 
ulator 

Uq,A)=exp{-{q/Af") (6) 
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Table 1. Skyrme force parameters. The column NLO list the 
parameters obtained from fits to the two-nucleon interaction 
Ref.[2], while the ones resulting from the fits to finite nuclei, 
Ref. [20) ■ are listed in column Skyrme. 



Parameter 


NLO 


Skyrme 


fo[MeV fm-"] 


-778 ... - 643 


-930 




0.17. . .0.41 


0.0... 0.40 


ti[MeY fm^] 


240 .. . 260 


77 . . . 130 


t2 [MeV fm^j 


-23 


-80. . .0 



for the Lippmann-Schwinger equation has been used. In 
nuclear matter, the states close to the Fermi surface dom- 
inate. Therefore, one can approximate the averaged con- 
tributions to the potential energy by replacing the mo- 
mentum in the regulator by the Fermi momentum and 
subsequently expanding the regulator, obtaining 

fn{q,A)v{q,q')Uq,A) = a + (3kp . (7) 

The term proportional to (3 can be interpreted as an ef- 
fective density-dependent interaction of the Skyrme type 
with a density dependence p", the power of the density 
dependence being a = 2n/3. A least squares fit to the 
nuclear ground state energies and radii for a large set of 
nuclei has been performed in Ref.fSOl as a function of the 
power a of the density dependence. The regulator em- 
ployed in Ref. [14] is defined by n = 2 and corresponds 
to an exponent a = 4/3 of the density dependence of 
the Skyrme interaction. The corresponding Skyrme pa- 
rameters to, xo,ti,t2 are displayed in Table [TJ The chiral 
two-nucleon interaction has contact terms proportional 
to Cs and Ct which translate into Skyrme parameters 
to = Cs — Ct and xq = 2Cs/{Cs — Ct)- Likewise, ti 
and t2 can be obtained from the NLO terms Ci and C2. 
Table [1] shows that the Skyrme parameters obtained from 
the two-nucleon interaction at order NLO agree in relative 
sign and are of the same order of magnitude as the ones 
deduced by fitting directly the bulk properties of finite 
nuclei. 

The present study demonstrates the applicability of 
the chiral two-nucleon interactions to nuclear matter and 
justifies the efforts required to investigate the next-to- 
next-to-leading order. 
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